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The one clean qubit model of quantum computation (DQC1) efficiently implements a computa¬ 
tional task that is not known to have a classical alternative. During the computation, there is never 
more than a small but finite amount of entanglement present, and it is typically vanishingly small in 
the system size. In this paper, we demonstrate that there is nothing unexpected hidden within the 
DQC1 model - Grover's Search, when acting on a mixed state, provably exhibits a speed-up over clas¬ 
sical with guarantees as to the presence of only vanishingly small amounts of quantum correlations 
(entanglement and quantum discord) - while arguing that this is not an artefact of the oracle-based 
construction. We also present some important refinements in the evaluation of how much entangle¬ 
ment may be present in DQC1, and how the typical entanglement of the system must be evaluated. 


I. INTRODUCTION 

Any computation whose quantum algorithm has a su¬ 
perior scaling of running time compared to its classical 
counterpart ought to pass through an intermediate state 
with non-trivial quantum properties, e.g. entanglement. 
Indeed, for pure-state computations, it has been shown 
that entanglement is necessary [1[. However, computa¬ 
tion involving mixed states is a far more subtle issue, 
with no firm resolution, although it is looking increas¬ 
ingly likely that quantum discord must be present in or¬ 
der for there to be an exponential speed-up (2}[3L 

The DQC1 model 0 can provide important insights. 
Expressed as a decision problem, this is defined as: 

Problem 1. Given an efficient classical description of a quan¬ 
tum computation U on n qubits, and a promise that either 
Tr(U) > l/poly(n) or Tr{U ) < —l/poly(n), determine 
which is the case with error probability e < |. 

There is no known classical algorithm which can effi¬ 
ciently solve this problem, while there is a quantum 
circuit that is remarkably simple (Fig. [Tb. In this cir¬ 
cuit, there is one special ('clean') qubit which is initially 
prepared in the state |+) = (|0) + \l))/\/2, and a set 
of n qubits are prepared in the maximally mixed state 
1/2". After applying controlled-!/ between the clean 
qubit (control) and the mixed qubits (target), we can 
estimate the probabilities that the clean qubit is in the 
state |+) (to find Re(Tr(t/))) or (|0) + i |1))/a/ 2 (to find 
Tm(Tr(7/))). There is a very obvious division in the sys¬ 
tem between the clean qubit and the mixed ones. As 
such, the authors of 0 examined the entanglement be¬ 
tween that bipartition of the system, and discovered it 
was 0. Entanglement seemed to not be necessary for 
mixed state quantum computation, although other mea¬ 
sures of non-classicality such as the quantum discord 
and measurement-induced disturbance have since been 
shown to be non-zero across that bipartition J5}|6j. 

Detecting entanglement in a multipartite state using 
bipartite divisions is, however, quite subtle. A clear 
demonstration of this is the distribution of entangle¬ 
ment using separable states |7] |8J; just because two of 



FIG. 1. The DQC1 model. There is no entanglement present 
in the bipartite split between the top (special) qubit and the 
rest. Estimation of the trace of U proceeds by performing 
measurements in the bases (|0) ± |l))/\/2 and (|0)±i|l))/ \/2. 

the three possible bipartitions of a 3-qubit system have 
no entanglement does not mean that the third biparti¬ 
tion has no entanglement. So it is with the DQC1 model 
- in EL it was demonstrated that by considering differ¬ 
ent bipartitions, there is indeed entanglement present. 
Two additional results are provided in 0: upper and 
lower bounds on the maximum amount of entangle¬ 
ment present across any bipartition and a numerical in¬ 
vestigation of the entanglement produced by a typical 
unitary (in the sense of selecting a U uniformly at ran¬ 
dom from the Haar measure). This investigation sug¬ 
gests that the entanglement is typically exponentially 
small in n. So, although there is some entanglement 
present, it is a vanishingly small amount on average. 

Even the presence of almost no entanglement in the 
computation might be seem surprising. However, in 
this paper, we draw parallels with a noisy version of 
Grover's Search ca in which we demonstrate the ex¬ 
istence of a quantum speed-up in the presence of van¬ 
ishingly small entanglement (such a result could have 
been proven from ||ll] but does not appear to have 
been), and other non-classicality measures such as the 
quantum discord lfl2l [T3l. Moreover, unlike the case of 
DQC1 where the speed-up over classical is only believed, 
Grover's search is subject to an oracle-based complex¬ 
ity classification with a proven gap over the classical 
case. Placed in this context, the power of DQC1 seems 
much less surprising. However, we only demand the 
existence of a quantum speed-up, and do not address 
the important conceptual transition between a polyno¬ 
mial and an exponential speed-up. We also improve the 
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results of J9), giving a tight upper bound on the entan¬ 
glement present in the DQC1 model for any U. 

A further criticism of the DQC1 model (or, similarly, 
the Grover search that we describe here) is that it as¬ 
sumes an implementation of the controlled-// gate. In 
practice, we have no such implementation, and need 
to decompose the action in terms of elementary gates. 
Even if the net implementation of controlled-// yields 
little or no entanglement, there is no a priori reason 
why all the intermediate steps should also be almost 
entanglement free across all Sipartitions. On the other 
hand, if the controlled-!/ is treated as an oracle, then we 
can consider other oracle based problems that are even 
more trivial in the entanglement sense - the Bernstein- 
Vazirani problem transforms a product state into a prod¬ 
uct state mm. The internal workings of the oracle are 
evidently important, and receive further discussion. 


A. Notation 

Consider a density matrix of n+1 qubits, p. Whenever 
we write p in the context of DQC1, we refer specifically 
to the output of the circuit for estimating the trace of U 
(Fig . [lj, which is a block matrix of the form 

l f t w\ 

P - 2 n+l [u t J ' 

We denote by p y , where y £ {0, l} n , the partial trans¬ 
pose of p taken over the non-clean qubits i : y, = 1, and 
never the clean qubit, without losing generality. Simi¬ 
larly, U y denotes the partial transpose of U, such that 

i_ /1 un 

Py 2 »+i t U, 1 J ' 

The Hamming weight of the string y is denoted by w y . 

Throughout this paper, the multiplicative negativity 
MM is chosen as the entanglement measure in order 
to maintain consistency with |9j. If A y = spec (p y ) is the 
spectrum of p y , then the multiplicative negativity of p is 

My = £ |A| 

AG A y 

for that specific bipartition. If no entanglement is 
present, M y = 1 for all y £ {0,1}". 


II. MAXIMUM ENTANGLEMENT IN DQC1 


Proof. Since it is not immediately clear that U y is normal, 
we use the singular value decomposition of U y , 

U y = RDV 1 

where R and V are unitary matrices and D is diagonal, 
with non-negative diagonal elements di . The eigenval¬ 
ues of p y are unchanged under the action of unitary op¬ 
erations, so apply the transformation of controlled-//' 
(controlled off the clean qubit), followed by a Pauli-A 
on the clean qubit, controlled-W, and finished by an¬ 
other Pauli-A' on the clean qubit. This yields 

_ 1 f 1 D\ 

Pv ^ 2n+1 ^ D 1 J • 

The eigenvalues of p y are therefore readily calculated to 
be (1 ± di)/2 n+1 for i = 1... 2 n . 

How, then, are we to pick {di} such that the value of 

My = l + 2 d 2 n + ! 

i:di> 1 

is maximised? We perform a constrained optimisation 
by noting from (9[ that 

2 n =Tr{UU^) = MU y Ul) = Y J dl 


Assume that t of the values d t > 1 (i = 1... /). Hence, 


M v < max 1- 

y ~ t,{di} 2 n 


1 

2 n 


I > 

i—1 


subject to i $ = 2". The optimal choice of the di is 
clear: di = 0 for i = t + 1... 2" and 


di — 



otherwise, which gives 

M v < max 1-—|- 

v ~ t 2 n 

This is maximised by / = 2"/4 (requiring n > 2 such 
that t is an integer), yielding 



We start by examining the entanglement properties of 
the DQC1 model, strengthening the assertions in (9). 

Theorem 1. For any unitary U, the output p of the DQC1 
computation ahvays satisfies 


for all possible bipartitions y £ {0,1}”. 


compared to the maximum possible value of ~ 2 n / 2 . □ 

The limit M y = | is readily saturated. For example, 
when n = 2 we can use 


U = 


0 0 0 1 \ 
0 10 0 1 
0 0 10 
1 0 0 0 / 
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since Uq\ has eigenvalues 2 , 0 , 0 , 0 as required for the op¬ 
timal construction. For larger n, a construction such as 
U 0 i ®{n- 2 ) would provide the requisite eigenvalues. 
There are then many ways of dressing this operator with 
unitaries that don't affect the partial transpose to dis¬ 
guise its structure slightly. For example, [9J uses a se¬ 
ries of controlled-not gates. Incidentally, this trivial ex¬ 
ample shows that there are many bipartitions for which 
M y = |: of the 2 n — 1 possible choices, all 2 n ~ 1 choices 
that have y\ 0 t /2 = 1 exhibit this value. 

This section reiterates the conclusion of 0 : there are 
bipartitions of the DQC1 model in which there is some 
entanglement present, and it can be present up to a fi¬ 
nite amount. Indeed, there can be many such biparti¬ 
tions. What we want to investigate in the rest of this 
paper is how surprising this result is - should the fact 
that we are guaranteed that no more than a small but 
finite amount of entanglement is present in a biparti¬ 
tion suggest to us that a quantum speed-up should be 
more difficult to realise? What about the fact that for 
most unitaries it seems that the entanglement is vanish¬ 
ingly small (exponentially small in the system size, N)1 
To that end, we are now going to compare to a depo¬ 
larised version of Grover's Search, for which we will 
show that even in situations where there is a provable 
(oracle-based) quantum speed-up, this can occur when 
we are guaranteed that there is never more than an ex¬ 
ponentially small amount of entanglement in any bipar¬ 
tition. This is a much stronger statement than has been 
made for DQC1. As such, using DQC1 to make conclu¬ 
sions about the existence of a computational speed-up in 
the presence of little or no entanglement seems obsolete. 


allowing \'4> ± ) to be expressed as 



Evidently, the state across this bipartition may be writ¬ 
ten as a state with no more than two Schmidt coeffi¬ 
cients. The maximum entanglement is when the two 
Schmidt coefficients are Ai = A 2 = 5 (w y = 1 as n — > 00 ), 
yielding a value of M y = 2 for this standard, pure state, 
version of Grover's search. There is only ever a finite 
amount of entanglement between any bipartition, much 
like DQC1, the only difference being the value. 

We will now show that there is a variant of this algo¬ 
rithm for which the entanglement can be made vanish¬ 
ingly small (no more than M y = 1 + 2 _T ^ 2 _e) for any 
e > 0). Consider using an initial state of 

p=p\+) <+r+d-p)i 

for any 0 < p < 1. A single run of the algorithm requires 
the standard quantum time, 0 ( 2 ra / 2 ), and succeeds with 
probability p = p + • By allowing L repetitions, there 

is a time advantage over classical (on average) if 2"/ 2 L < 
2 n_1 , and it succeeds with probability 1 — (1 — p) L ss 

pL for small p. As such, for any finite e, p ~ 2 _ "^ 2 _e i 
provides an advantage over classical search algorithms 
while we evidently anticipate that for smaller p, there is 
less entanglement present. 


III. COMPARISON WITH GROVER'S SEARCH 


A. Entanglement 


A standard formulation of Grover's Search algorithm 
HTDl starts with an initial state |+)® n , trying to evolve 
to some (unknown) target state |aro)- The evolution is 
restricted to a subspace spanned by \xq) and 


VO 


1 

sj2 n - 1 


I s )’ 

X^xo 


such that the populated states are of the form 

I’T) = cost* |xo) + sin# \ip ± ) 


We turn to calculating the entanglement present in the 
depolarised algorithm. Progress through the algorithm 
can be specified by a parameter 9, such that the state is 

p(0)=p|^)) W)| + (l-p)^ 

Assume the Schmidt coefficients for the state T) are 
cos 2 (j> and sin 2 <j> with respect to bipartition y. In the 
large n limit, <f> can take on any value ( 0 , 7 r/ 2 ]. p(9 ) is 
unitarily equivalent (where the unitaries are local with 
respect to the bipartition, and hence do not change the 
eigenvalues under the partial transposition) to a state 


where 9 G [0, 7 r/ 2 ], For finite n, only integer multiples of 
sin 1 (2 _ ") are allowed. If we select a particular bipar¬ 
tition, then |cco) = | x y ) \x y ) where, here, the subscript y 
merely refers to the set of qubits on which the state is 
defined. Similarly, we can define 


„ ._ —_ y 

v 1 W2 w y - 1 ^ 

xE{0,l } w y :x^x y 


\K) = 


( p cos 2 <j> + 0 0 p cos <f> sin </> ^ 

0 ' y 0 0 

0 0 y 0 

p cos (j) sin cj) 0 0 p sin 2 <f> + 

V ^l2"-4 ) 

( 1 ) 

where the upper left 4x4 describes the space spanned 

by \x y )\ Xy), \Xy)\l/)±), \i>y)\Xy) &Ud \^y)\^)- ^ 
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can readily take the partial transpose of this, and cal¬ 
culate the eigenvalues: (repeated 2" - 4 times), 

|(1 ± cos(2</>)) + and ± § sm(24>). Hence, 

M v = 1 ^ + psin(2</>). 

This is maximised at (f> = n /4, indicating that for all bi¬ 
partitions, and at all points during the computation, 

m < i+p- r- 
2" -1 

When p is exponentially small, the entanglement is al¬ 
ways exponentially small and there is still a provable 
(oracle-based) speed-up in searching over the classical 
case. This compares favourably to DQC1 wherein typ¬ 
ical unitaries have been numerically shown to produce 
exponentially small amounts of entanglement (9)/ while 
having a believed speed-up over classical. 


B. Quantum Discord 


Perhaps it is not entanglement that needs to be 
present in a mixed state quantum computation. Instead, 
other measures of non-classicality, such as the quantum 
discord have arisen. If exponentially small entangle¬ 
ment might be considered surprising, perhaps it is the 
case that there is finite quantum discord? We will now 
show that the discord is also vanishingly small. A pre¬ 
cise definition Ifl2l [13 ] of the discord is unnecessary. It 
suffices to know that it is a non-negative quantity which 
is always 0 in the classical limit and, for pure states, re¬ 
duces to the entanglement entropy. 

We take p(ff) as for the previous calculation. 


/ p cos 1 2 (j) + 0 0 p cos </> sin (j> 

0 ^0 0 
0 0 ^ 0 

p cos (j) sin <j> 0 Op sin 2 </> + 

V 

and define a second state <j = 

/ p cos 2 (j) + 0 0 0 

0 ^ 0 0 

0 0 ^ 0 

0 0 0 p sin 2 <j> + 

V 


\ 

^l 2 n _ 4 / 

\ 


By virtue of being diagonal with respect to a separable 
basis, the state a has 0 discord m, and yet is very close 
in terms of trace distance to p: 


d = Tr|p — a | = psin(2</>). 

The quantum discord is continuous Ifl9l [20], meaning 
that we can bound the amount of discord present in p(9): 

D(p) < 8 d{n — 1) + 4 h{d) 


where h(x ) = — x log 2 x—(l — x) log 2 (l — x) is the binary 

entropy. This scales as 0{n2~ n ^~ e ^), and is therefore 
also vanishingly small across all bipartitions simultane¬ 
ously, and for all steps of the algorithm. 


C. Decomposing the Oracle 

An oracle-based problem (whether this is explicit, as 
for Grover's search, or implicit, requiring the implemen¬ 
tation of controlled-t/ in DQC1) is useful for determin¬ 
ing the computational complexity of a problem (with 
respect to that oracle), often enabling lower bounds as 
well as upper bounds. However, it has the potential to 
mask a lot of the entanglement properties. If one has 
to implement an oracle by using a collection of smaller, 
more manageable, operations, we need to be sure that 
each of those individual operations does not leave an 
intermediate state that has large amounts of entangle¬ 
ment. After all, there are many quantum computa¬ 
tions that start and end in separable states, and it is 
only the intermediate products that are entangled. The 
Bernstein-Vazirani algorithm is one case that hides all 
the entanglement in such a way 114, 15|. Typical DQC1 
computations probably also hide their entanglement - 
imagine decomposing U in terms of one-qubit unitaries 
and controlled-nots, and implementing the controlled-t/ 
by consecutively applying the appropriate controlled- 
one-qubit and Toffoli gate^J The first Toffoli gate is 
likely to introduce a finite amount of entanglement. 

We claim that there need be no corresponding diffi¬ 
culty when using Grover's search. The Grover iterator 
may be written as 


H ®np Q H® n P X0 , 

where H is the Hadamard gate, and P x is an ri-qubit 
controlled-phase gate which adds a phase of 7r only to 
the state \x). Let us take each step in turn. If p m is the 
state before the action of the iterator, and p out is the state 
after the application, then our results so far show that 
both have exponenitally small entanglement and dis¬ 
cord. Now, observe that after the action of P Xo , the form 
of pin in Eq. ({Th is only changed by adding a negative 
sign on the off-diagonal elements. Both the entangle¬ 
ment and discord are unchanged. Application of a set 
of local unitaries (H® N ) also cannot change these val¬ 
ues. What about the application of PqI Instead, we ob¬ 
serve that after the action by P 0 , this is the same as state 
Pout with H® n applied to it, which must have the same 
entanglement properties as p out . 

Should we decompose the operations any further, 
perhaps by writing P XQ in terms of a universal set of 


1 This is certainly far from the only decomposition that one could 

make, and this is merely an illustrative argument. 
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two-qubit gates? We suggest this is unnecessary as it 
appears that in a wide variety of experimental imple¬ 
mentations, the gate P Xo can be implemented directly 
62 TH 23 I . So, we are justified in our claim that every el¬ 
ementary step of the algorithm leaves us in a state of 
almost no entanglement and almost no discord. 

IV. CONCLUSIONS 

One of the aspects of the DQC1 model of compu¬ 
tation that originally generated much interest was the 
suggestion that it achieved its computational speed-up 
without entanglement. In fact, it does use entanglement 
0, and that entanglement can be up to a finite amount 
across many different bipartitions of the system. Even 
on those bipartitions for which there is no entanglement, 
it has been shown |5j instead that other measures of 
non-classicality are present in finite amounts and could 
be a resource for the computational speed-up that is be¬ 
lieved to be present in DQC1. 

In this paper, we have contrasted this with Grover's 
search when acting on an initial state that has a large 
admixture of the maximally mixed state. The degree 
of admixture for which a computational speed-up is 
still possible (without any attempt at error correction) 
was readily derived and, as such, it was shown that 
over every bipartition both the entanglement and quan¬ 
tum discord are vanishingly small. Moreover, we know 
that the computational speed-up is present for Grover's 
search, since we know the minimum running time for 
the classical algorithm (using a certain oracle). We con¬ 
clude that DQC1 is nothing special - Grover's search ex¬ 
hibits stronger properties in every way (guaranteed c.f. 
assumed speed-up, guaranteed exponentially little en¬ 


tanglement c.f. little entanglement on average, exponen¬ 
tially little discord c.f. finite discord), except for one fea¬ 
ture - DQC1 is believed to have an exponential speed¬ 
up while Grover's search is only polynomial. 

The presence of quantum discord in the DQC1 model 
has been taken by a number of authors as "the first 
real evidence that mixed-state quantum computation 
can have an advantage over classical computation even 
when entanglement is absent" | |24| . choosing to view the 
almost-no entanglement as essentially equivalent to no 
entanglement. However, here, we have almost no entan¬ 
glement, and almost no quantum discord (and, indeed, 
any measure of quantum resource that is continuous). 
We are left to conclude that there is a marked difference 
between 'almost no' and 'no' entanglement/discord. 
Exactly what is required for a computation to gain speed 
over a classical one remains unclear, but it seems that 
only small amounts of whatever resource suffice. How¬ 
ever, this leaves open the possibility that larger amounts 
of a resource, such as quantum discord, are necessary for 
an exponential vs. polynomial speed-up. 

On a final note, we wish to raise an important issue 
with the calculation of the typical entanglement 0 or 
discord |5] across a bipartition in the DQC1 computa¬ 
tion. The decision problem for DQC1 was defined by ex¬ 
plicitly bounding the trace away from 0. However, typ¬ 
ical unitaries have an expected trace of 0, and the trace 
is closely centred on the 0 value. As such, the decision 
problem explicitly asserts that the unitaries involved are 
not typical, and further studies are required to clarify 
the impact of this, although the studies of special cases 
that are performed in the Appendix suggest that the im¬ 
pact is not significant. 

Acknowledgements: We thank L.C. Kwek for introduc¬ 
tory conversations. 
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Appendix A: Typical Entanglement 

In this Appendix, we illustrate some important fea¬ 
tures of the calculation of typical entanglement in the 
DQC1 model. Rather than permitting an arbitrary Haar- 
random unitary, we concentrate on a subset of unitaries, 
based on graph states, that facilitate simpler calcula¬ 
tions. We define the random unitary, U, in the following 
way: consider n qubits and select a random diagonal 
matrix, meaning that each diagonal element is selected 
at random to be e >f>: ' for all x G {0,1}™. We then conju¬ 
gate this with Hadamard gates on each qubit (i.e. both 
before and after the diagonal matrix) before conjugat¬ 
ing by controlled phase gates between all pairs of qubits 
(i, i + 1). These conjugations are the unitaries required 
for preparing the one-dimensional cluster state. 

The reason for concentrating on graph states is that 
there is a well established formalism for calculating the 
partial transpose Il8l l25ll26l . If the eigenvalues of U are 
written as a vector 


i©> = E eie * !*>’ 

X 

then the eigenvalues of E/01010101...01 are given by the 
vector R\Q) where 

R = it®” ( E XiXi+1 I*) {x\ ] H® n . 

\xe{o,i}« ) 

Lemma 1 . If n is even, then every matrix element of R has 
| (y\R \z) | = In every rozv, there are ^(2 n + 2 n / 2 ) 

positive entries, and |(2" — 2"/ 2 ) negative entries. 

Proof. The matrix elements have the form 

{y\R\z) = ^ E (-ir <z+v) (-^^ XiXi+1 - 

xe{o,i} n 

First, observe that (y\ R\z) = (000 ... 0| R\z ® y), so it 
suffices to concentrate solely on the first row. The ele¬ 
ments of all other rows are just permutations of the first. 
We will denote the elements of this first row by R z . 
Notice that 


E R, = 1. 

ze{o,i}" 

This means that if it is true that all 2" elements have 
\R Z | = 1 /2"/ 2 , then since they are all real, it must be that 
there are \{2 n + 2 ra / 2 ) positive entries, and |( 2 " — 2 n / 2 ) 
negative entries. 

We will prove the first half of the lemma by induction, 
using the base case of n = 2: 

f?(00) = \ f?(01) = i P(10) = R{ 11) = 5 . 

Assume that |i?( n_ 2 )(z)| = 2 „ / 1 2 _ 1 . Now con¬ 
sider R( n fz\\z n -iz n ). For each term x in the sum 


0.25 
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FIG. 2. Comparison of theoretical results (lines), Eq. 
and average for 10000 samples (points) for both ID (upper, 
blue) and 2D (lower, purple) random walks. 

for R^ n ~ 2 \z), there are 4 terms in R^ (z\\z n -\z n ): 
x||x„_i x n . We can consider each of them separately, and 
how they affect the additional phases 

Z n X n-\-Zn — lX rl -l-\-X rl X rl — l-\-X rl — iX rl — 2 

One of the phases can be moved into an effective z: 

R^(z\\z n _ lZn )=^ e 

X n ,X n — l 

R { - n ~ 2 \z®{ 00...0a; n _i)). 

On performing the sum over x n , we either have 
iR^^fz) (if = 0 ) or ±R( n - 2 Xz © ( 00 ... 01 )). 
By assumption both have the same absolute value, 

l/2”/ 2 . □ 

Our aim is to evaluate the multiplicative negativity: 

M y = 1 + frv ( A-1 )- 

Aefl|0):A>l 



We assume that each element of R |0) is independent 
(this may not be strictly true, but is a good approxima¬ 
tion). This reduces M y to M y = 1 + E where E is the ex¬ 
pected value of max(A — 1,0). To evaluate this, we need 
the probability distribution for A. As a consequence of 
the Lemma, we have 



f 

i (2 "+2"/ 2 ) 2 n 

\ 

P(A > 1) = P 

\ 

E E 

i - 1 i=i( 2 "+ 2 "/ 2 ) + l 

> 2 n / 2 


However, if the probability distribution of 0, is unaf¬ 
fected by a 7 r shift, this is entirely equivalent to 


P 





We consider two cases. Firstly, each 0i is a random 
choice, 0i € {0, 7 r}. In this case, (E=i el8i i s simply 














7 


* • 

• • 

—•—*-*-*- - -Number of Qubits 

s in is in 


FIG. 3. Typical entanglement for a graph diagonal unitary of 
fixed trace Tr([/)/2" = §, |, | for shades black through to 

light grey respectively. The larger the trace, the less entan¬ 
glement is present. 


the expected distance of a random walk in ID. Secondly, 


9i £ [0, 27t) leads to the interpretation of 



as a 


random walk in 2Erl In either case, Rayleigh's solution 
for the probability distribution at large n is 


P 




< x + Sx 


ID 

§f e-^'^dx 2D 


in order to calculate 


E = 





< x + Sx 


X 

2 n /2 



which yields 


E= l 1D 

1 ^erfc(l) 2D 


(Al) 


Fig. [2] provides numerical confirmation of this calcula¬ 
tion. vVe conclude that in both cases there is a constant 
amount of entanglement. This is in contrast with the nu¬ 
merical results in 0 for general random unitaries. 

However, there is one further subtlety that it is im¬ 
portant to raise. Recall that for the DQC1 model we are 
given the promise that the trace of U is bounded away 
from 0. Meanwhile, the probability distribution for the 
trace of a typical unitary is strongly centred on the value 
0 - in effect, we are post-selecting on highly atypical 
unitaries so typicality arguments might not reveal ev¬ 
erything. However, it turns out that this only serves to 
reduce the amount of entanglement present. For a sim¬ 
ple argument, consider the case where Tr(f7) = 2 n . In 
this case, we know that U = 1, and that M y = 1, less 
than the above-calculated values. More generally, we 
have numerically examined the random unitary model 
described above (restricting the eigenvalues to ±1), fix¬ 
ing the trace to be different values of Tr(f7), see Fig. [ 3 ] 
The amount of entanglement present decreases as the 
value of the trace is increased. 


2 We define this to be the walk such that at each step, a step length of 
1 is taken in a random direction in the plane. Some authors choose 


to define it as a walk on a square lattice. 


















